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Abstract 

The QCD evolution of the pion distribution amplitude (DA) ^(x, Q 2 ) is computed for several 
commonly used models. Our analysis includes the nonperturbative form predicted by light-front 
holographic QCD, thus combining the nonperturbative bound state dynamics of the pion with 
the perturbative ERBL evolution of the pion DA. We calculate the meson-photon transition 
form factors for the 7r°, 77 and rf using the hard-scattering formalism. We point out that a 
widely-used approximation of replacing <j) (x, (1 — x)Q) with Q) in the calculations will un- 
justifiably reduce the predictions for the meson-photon transition form factors. It is found that 
the four models of the pion DA discussed give very different predictions for the Q 2 dependence 
of the meson-photon transition form factors in the region of Q 2 > 30 GeV 2 . More accurate 
measurements of these transition form factors at the large Q 2 region will be able to distinguish 
different models of the pion DA. The rapid growth of the large Q 2 data for the pion-photon 
transition form factor reported by the BABAR Collaboration is difficult to explain within the 
current framework of QCD. If the BABAR data for the meson-photon transition form factor 
is confirmed, it could indicate physics beyond-the-standard model, such as a weakly-coupled 
elementary C = + axial vector or pseudoscalar z° in the few GeV domain, an elementary field 
which would provide the coupling 7*7 — >► z° — >> tv° at leading twist. Our analysis thus indicates 
the importance of additional measurements of the pion-photon transition form factor at large 

Q 2 . 

PACS numbers: 13.40.Gp, 12.38.Bx, 14.40.Be, 14.40.Df, 13.60.Le 
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I. INTRODUCTION 



The BABAR Collaboration has reported measurement of the photon to pseudoscalar- 
meson transition form factors from the 7*7 M process for the 7r° [I], r\ c [2], rj, and 
rf [SI Hj. The momentum transfer Q 2 range covered by the BABAR experiments is much 
larger than the range studied by the CELLO [5j and CLEO [6j collaborations. More 
significantly, the BABAR data for the 7T°-7 transition form factor exhibit a rapid growth 
for Q 2 > 15 GeV 2 which is unexpected from QCD calculations, whereas the data for 
the other transition form factors agree well with previous measurements and theoretical 
calculations. 

QCD computations for exclusive processes are considerably more subtle than inclusive 
processes since one deals with hadron dynamics at the amplitude level. The foundation 
for calculating exclusive processes at high momentum transfer in QCD was laid down 
almost 30 years ago [7H9]. It was shown in [7j that the pion electromagnetic form factor 
and transition form factor (TFF), the simplest exclusive processes involving the strong 
interaction, can be calculated as a convolution of a perturbatively-calculable hard scatter- 
ing amplitude (HSA), and the gauge-invariant meson distribution amplitude (DA) which 
incorporates the nonperturbative dynamics of the QCD bound-state. The distribution 
amplitude, Q) is the qq light-front wavefunction (LFWF) ^(x,k_L), the eigenstate 
of the QCD light-front Hamiltonian in light-cone gauge, integrated over transverse mo- 
menta < Q 2 . Here x = k + /P + = (k° + k z )/(P° + P z ) is the light-front momentum 
fraction of the quark. The DA has the physical interpretation as the amplitude to find 
constituents with longitudinal light-front momentum x and 1 — x in the pion which are 
non-collinear up to the scale Q. The form of the DA can be confronted with the results of 
various processes sensitive to the form of the DA and calculated using non-perturbative 
methods [ID]. There are also important constraints from the lowest moments of the pion 
DA obtained from lattice gauge theory [IIJ [12] . 

The evolution of the pion DA is governed by the Efremov-Radyushkin-Brodsky-Lepage 
(ERBL) equation [7-9 J . The form of the pion DA: cf) asy (x) = V3f n x(l-x) and the result- 
ing predictions for elastic and transition form factors at the asymptotic limit Q 2 — >■ oo 
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can be predicted from first principles [7J . The results are independent of the input form 
of the distribution amplitude at finite Q 2 . However, the prediction for the elastic form 
factors using the 'asymptotic' form for the pion DA at finite Q 2 range are not successful 
when compared with available experimental data. This has led to many theoretical inves- 
tigations of the shape of the pion DA at low Q 2 which reflect the bound state dynamics. 
Some models which are vastly different from the asymptotic form have been suggested; 
however, these forms lack physical motivation and contradict the lattice constraints. This 
is manifested by the suggestion of a 'flat' form [T3l - fT5] for the pion DA in order to explain 
the recent BABAR measurements pQ for the pion TFF. 

The effects associated with the transverse momentum degree of freedom have been 
analyzed in Refs. [T6HT8] . It was shown in [16] that the transverse momentum dependence 
in both the HSA and the LFWF needs to be considered in order to make predictions for 
the pion TFF for Q 2 of the order of a few GeV 2 . 

The pion form factor has been calculated using the asymptotic DA and Chernyak- 
Zhitnitsky (CZ) form [19J at next-to-leading order (NLO) [201423] . using the standard 
hard-scattering approach when the k^-dependence in the HSA is ignored. The next-to- 
next-to-leading order corrections to the hard-scattering amplitude were calculated in [24] 
using the conformal operator product expansion. The form factor has also been studied 
[T71 [T8] using the modified hard scattering approach in which the k^-dependence is 
considered together with gluon radiative corrections. In these calculations the evolution 
effects were often shown together with high order corrections. However, due to the 
limitation on the form used for the pion DA, the effects from evolution have not been 
fully explored. There are many other theoretical studies of the pion-photon transition 
form factors (see for example [25HIU] ). 

In this paper, we reexamine the relation between the light-front wavefunction and the 
distribution amplitude and calculate the meson-photon transition form factors for the 
7T°, 7] and ?/. Various forms of the meson distribution amplitude and their evolution 
are studied in Section II. Our analysis integrates the nonperturbative bound state dy- 
namics of the pion predicted by light-front holographic QCD with the perturbative QCD 
ERBL evolution of the pion distribution amplitude, thus extending the applicability of 
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AdS/QCD results to large Q 2 . The pion-photon transition form factors for the real and 
virtual photons are calculated in Section III. The ?7-photon and ry'-photon transition form 
factors are studied in Section IV. Some conclusions are given in the last section. 



II. PION LIGHT-FRONT WAVEFUNCTION AND DISTRIBUTION AMPLI- 
TUDE 



The pion distribution amplitude in the light-front formalism [Zj is the integral of the 
valence qq light-front wavefunction (LFWF) in light-cone gauge A + = 



(x,Q)= / 
Jo 



d 2 k ± 
16vr 3 



(1) 



The pion DA can also be defined in terms of the matrix element of the axial isospin 
current between a physical pion and the vacuum state [H] 

(Q) 



<f>(x,Q) 



dz 
~2tt~ 



J(2x-l)z-/2 / q 



7T 



z+=z ± =0;pt=0 



(2) 



where 



fl = exp |ig J dsA + (zs)z 



(3) 



is a path-ordered factor making Q) gauge invariant. The pion DA satisfies the 
normalization condition derived from considering the decay process tt — >■ \±v [Nq — 3) 



I 



1 U 

da; <p(x, n) 



(4) 



/o 2^3' 
where = 92.4 MeV is the pion decay constant and /i is an arbitrary scale. 

By definition the bound state LFWF ^^(rr, kj_) has important support only when 
the virtual states are near the energy shell, i.e. 



x(l — x) 



< /4> 



(5) 



where \±p can be viewed as the factorization scale. Thus a 'cut-off' on the transverse 
momentum is implied in the definition for the soft component of the LFWF: '0^/ t 7r (^, kj_). 
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A natural way to implement this cut-off is to require the LFWF to decrease quickly for 
large k^, for example, via an exponential function as first suggested in the model discussed 
in [42]. One can write a parameterization form for the LFWF as in [43j 

^(x,k ± ) = #r)£(x,k ± ) 
where k is the gap parameter, and 



and the function E satisfies 1 , 

r0 ° d 2 k i 



-£(x,k ± ) = l. (8) 



16tt 3 ' 

A common practice used in the literature in determining the parameter k is calcu- 
lating the non-perturbative properties of the pion and comparing with the experimental 
measurements of these quantities. However, this process only allows one to constrain 
k in a relative large range due to the uncertainty of the experimental measurements. 
For example, the root of the mean square transverse momentum of the valence quarks, 
defined as 



2\ 



where P q q is the probability of the valence Fock state of the pion 



qq 

-1 poo j2 



pi poo 72k. 

P« = I dx j( j^\^%(x,k ± )\ , (10) 

is estimated to be in the range of 300 ~ 500 MeV from experimental measurement on 
the charge radius of the pion. Thus k is not well determined by Eq. ([9]). 

Brodsky, Huang, and Lepage [42] obtained a constraint for the soft LFWF at = 0, 
^^/tt( x ' k-L = 0), by studying the decay of 7r° — >> 77. However, we note that the decay 



1 Strictly speaking a cut-off of |k^_| max ~ x(l — x)n 2 F is still in place for the soft wave function given by 
Eq. (|6j). However, calculations are not sensitive to this cut-off due to the nature of rapid decreasing of 
the wave function. Thus it is commonly expressed in the literature that |k^_ | max = ji 2 Fl or |k^_ | max = 00. 
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7T° —> 77 is a long-distance process for which the higher Fock states should make substan- 
tial contributions as well, since there are extra interactions with the quark propagator 
between the two photons which vanish at high Q 2 in the light-cone gauge. Therefore k 
cannot be well determined from the constraints imposed by the decay process. 

From these considerations we will treat k in Eq. ^ as a phenomenological parameter 
which is allowed to change in a certain range. It is equivalent to treat the probability of 
the valence Fock state, P qq) or the root of the mean square transverse momentum of the 
valence quarks, \J (kj_), as a parameter. 

The LFWF i/j qq / 7r (x,\£±) in Eq. ([!]) contains all the non-perturbative information of 
the pion. There are also perturbative corrections that behave as a s (k^)/k^ for large k^, 
coming from the fall-off of the LFWF ^(x, kj_) due to hard gluon radiation [3, [42]. Both 
soft and hard regimes are important to compute the pion transition form factor for all 
values of Q 2 . 



A. Soft Evolution of the Pion Distribution Amplitude 



Substituting Eq. @ into Eq. Q one obtains 



Q) = </>(x) 



exp 



Q 2 



ii) 



2k 2 x{\ — x) 

where (j)(x) is given by Eq. Q. Eq. (11) gives a factorization model for the Q 2 dependence 



of the distribution amplitude in the soft domain. The soft Q dependence in Eq. (11) can 



be safely ignored for Q > 1 GeV for the typical values of k ~ 0.5 — 1.0 GeV. In the regime 
of Q > 1 GeV one needs to consider the hard gluon exchanges that provide additional 
logarithmic Q 2 dependence in Q), as given by the ERBL evolution equation discussed 
below. 

Many efforts have been made in determining the pion DA at a low momentum transfer 
scale /ig ~ 0.5 — 1 GeV. Most of these studies concentrate on the determination of the 
first few terms in the solution of the evolution equation for the pion DA discussed in 
the next section. However the pion DA at a low scale could differ significantly from its 
asymptotic form due to the slow convergence of the evolved DA. Using only a few terms 
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of the full solution will put a strong limitation on the studies. The following forms have 
been suggested. 



(a) The asymptotic form UU PB] 

< f> as y(x) = V3f„x(l-x). (12) 



(b) The AdS/QCD form 03 M 



AdS (x) = 4 /7rV ^T3^ ). (is) 

V37T 



(c) The Chernyak-Zhitnitsky [19J form 

cz (x) = 5^3/^(1 -x)(l-2x) 2 



l + |cf 3 >(l-2s) 



(14) 



(d) The 'flat' form [[5 



/ at (x) = ^= [TV + 6(1 - - x)} . (15) 

The DA model (b) follows from the precise mapping of string amplitudes in Anti-de 
Sitter (AdS) space to the light-front wavefunctions of hadrons in physical space-time 
using holographic methods [451448] . However, an extended AdS model with a Chern 
Simons action maps to the asymptotic DA form x(l — x) [49j rather than the AdS form 
y/x(l — x). A discussion of the pion form factor is discussed in the framework of light- 
front holographic mapping in a forthcoming paper [50j . Model (c) was suggested on the 
basis of a calculation using QCD sum rules and model (d) was advocated in explaining 
the recent BABAR data for the pion TFF pp. The end-point non- vanishing models, similar 
to model (d), were also obtained [5TH53] for the pion and photon DAs using chiral quark 
models and Regge models before the BABAR results were reported. Normally one expects 
that the light-front wavefunction of a composite hadron to vanish at the x — 0, 1 end- 
points to ensure a finite expectation value of the kinetic energy operator. A set of pion 
DAs (termed the BMS models) including only the first two terms in the general solution of 



the ERBL (see Eq. 16) below) were proposed [25] by comparing the light-cone sum rules 
calculations for the pion TFF with the CELLO and CLEO data. Theoretical calculations 
using transverse lattice gauge theory with discrete light cone quantization [54] and chiral 
quark models [55] generally suggest that that the pion DA is considerably broader than 
the asymptotic form. 

0.6 
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FIG. 1: The four commonly used models for the pion distribution amplitude. The curves from 
bottom to top at x — 0.5 are for the CZ, 'flat', AdS/QCD, and asymptotic forms, respectively. 



The four models are shown in Fig. [TJ Model (d) is not actually only flat over the 
whole range of x - it is end-point enhanced. Models (c) and (d) have very different shape 
from (a) and (b). The zero- value of the CZ form in the middle point (x = 0.5), where 
the pion momentum is shared equally between the quark and the antiquark, and the 
enhancement of the 'flat' form in the end-points (x = 0, 1), where the pion momentum 
is mostly carried by the quark or the antiquark, are hard to understand in terms of the 
bound state dynamics of the pion. The zero-value of the CZ form in the middle point 
also disagrees with the estimation using the QCD sum rule method reported in [56] . 
<f>(x = 0.5) = (0.17 ±0.03)/*. 

Using the models of the pion DA discussed above we can construct the corresponding 
LFWF from Eq. The LFWF constructed with the 'flat' form (model (c)) for the 
pion DA is non-normalizable since the probability of finding the valence Fock state in 
the pion (Eq. (p~0|) ) becomes infinity [55] . We list the values for the gap parameter k and 



the root of the mean square transverse momentum of the valence quarks \J (kj_) for the 
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TABLE I: Properties of the soft light-front wavefunction corresponding to various models of 
the pion DA. 



p - 


DA 


k (GeV) 


V< k i> ( GeV ) 






0.826 


0.370 


0.25 




0.859 


0.350 




</> cz (x) 


1.210 


0.403 






0.584 


0.261 


0.50 


</. AdS (x) 


0.607 


0.248 




cf> cz (x) 


0.855 


0.285 






0.462 


0.207 


0.80 


4> AdS (x) 


0.480 


0.196 




</> cz (x) 


0.676 


0.225 



three choices of the probability P q q = 0.25, 0.50 and 0.80 in Table [IJ For the 'flat' model 



Eq. (10) is divergent so we adopt k 2 = 0.530 GeV 2 [Hj and N — 1.3 [T5] . which were 
chosen to explain the BABAR data. 



B. Hard Evolution of the Pion Distribution Amplitude 



The evolution of the pion DA at large Q is governed by the ERBL equation. The so- 
lution to the ERBL equation can be expressed [3 [8] in terms of Gegenbauer polynomials, 



<f>(x,Q)=x{l-x) Yl a n (Q)C 3 J 2 (2x-l), 

71=0,2,4,- 



(16) 



where 



a>n(Q) 



(17) 



[a s (Q 2 )_ 

at leading order. The coefficients {a n (/x )} are the coefficients in the Gegenbauer expan- 
sion of the DA at the initial scale /xq, 



(x, fjio) = x(l -x) ^2 an{^)Cl /2 {2x - 1), 

n=0,2,4,- 



(18) 
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and follow from the orthonormality of the Gegenbauer polynomials 

anM = i ^S^Tn [ 1 dx4 > (x,^)C 3 J 2 (l-2x). (19) 
(n + 2)(ra+ 1) J 

The QCD coupling constant a s (Q 2 ) is taken to have the leading-order form 

A)ln (<57 A qcd) 

where Aqcd is the QCD scale parameter and (3q is the QCD beta function one-loop 



coefficient flo = 11 — {^rif. The anomalous dimensions 7 n appearing in Eq. (17) 

n+l 



_ 4 

In- 3 



3 + 



" T1 1 

4 E7 



n + l)(n + 2) 



(21) 



are the eigenvalues of the evolution kernel |8]. The coefficient a (/x ) = V^/tt for 
any model of the pion DA, since the pion DA should satisfy the normalization condition 
Eq. Q with Cl ,2 (z) = 1, and - x)Cn /2 (l - 2x) = for n > 2. 

The coefficients a n (/x ) are computed at the initial scale /x = 1 GeV (where the 
effects of hard gluons is negligible and the scale dependence of the soft evolution is not 
important). Thus we choose the initial condition 0(x,/x o — 1 GeV) ~ with (f>{x) 

given by Eq. 0. At leading order the asymptotic form does not evolve since all the 
expansion coefficients {a n (/x )}, but a (/x ) = vanish for model (a). The coefficients 

{a n (/xo)} for model (c) (the CZ form) are very simple since the model essentially includes 
only the first two terms in the Gegenbauer polynomials, a (/io) = V^fn and a 2 (/x ) = 
2/\/3 f n . For model (b) (the AdS form) we include the first 50 terms (i.e. up to n = 100) 



in Eqs. ( |T6| ) and (18) in our calculation. The first 10 values of a n (/xo) for the AdS model 
and the nonzero coefficients for the asymptotic and CZ models are listed in Table [TTJ It 
was found that for the AdS model the calculation with 51 terms only brings a few percent 
corrections to the calculation with 21 terms over a large range of x. 

It is problematic to expand the 'flat' DA in term of the Gegenbauer polynomials at the 



initial scale /x , since the expansion Eq. (18) converges if, and only if, 0(x,/x o ) vanishes 
at end-points [3, [57] . We will not try to apply the ERBL equation to the 'flat' DA, but 
just make the note that if one applied the ERBL equation to the 'flat' DA, one would 
enforce the suppression at the end-points as soon as the evolution starts. 
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TABLE II: The coefficients a n (/io) for the asymptotic, AdS and CZ models for the pion DA. 



n 





2 


4 


6 


10 


12 


14 


16 


18 


20 


On(Mo)/ao(Mo) 


AdS 


1 


0.1461 


0.0573 


0.0305 


0.0189 


0.0129 


0.0094 


0.0071 


0.0056 


0.0045 


CZ 


1 


2/3 


for n > 4 


asy 


1 


for n > 2 



The first term in Eq. (16) represents the asymptotic form of the pion DA and the 
asymptotic form does not evolve with Q 2 . The other distribution amplitudes have Gegen- 
bauer polynomial components with nonzero anomalous dimensions which drive their con- 
tributions to zero for large values of Q. One can start with any distribution amplitude 
0(x,/io) at any finite scale and expand it as x(l — x) times Gegenbauer polynomials. 
Only its projection on the lowest Gegenbauer polynomial with zero anomalous moment 
survives. This is illustrated in Fig.[2]for the first few expansion coefficients a n (Q 2 ) for the 
AdS distribution amplitude. The evolution effects of the DA at leading order are shown 




50 100 



G 2 (GeV 2 ) 



FIG. 2: Evolution of the expansion coefficients a n (Q 2 ) for the AdS distribution amplitude. The 
curves from top to bottom are for n = 2, 4, 6 and 8 respectively. 

in Figs. [3] and [4] for the AdS model and CZ model for the pion DA respectively. In our 
numerical calculations we used fio = I GeV and Aqcd = 225 MeV. Performing evolution 
at NLO modifies the results slightly. It can be seen that evolution effects change the 
shape of the CZ form significantly, while the effect on the AdS form is not as dramatic. 
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In the asymptotic Q 2 — >> oc limit the asymptotic DA is recovered. 




x 



FIG. 3: Evolution effects shown for the AdS model for the pion DA. The curves from bottom 
to top at x — 0.5 are for Q 2 — 1, 10, 100 and 1000 GeV 2 , and the asymptotic DA, respectively. 

C. Moments of the Pion Distribution Amplitude 

Important constraints for the form of the distribution amplitudes also follow from 
QCD lattice computations. The latest results for the second moment of the pion DA, 




0.0 0.2 0.4 0.6 0.8 1.0 



x 

FIG. 4: Similar as in Fig. [3] but for the CZ model for the pion DA. 
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where f = 1 - 2x, are (£ 2 ) M 2 =4GeV 2 = 0.269 ± 0.039 [II] and 0.28 ± 0.03 [E]. The second 
moments calculated at the initial scale /ig for the four models of the pion DA described 
above are 0.20 (asymptotic), 0.25 (AdS/QCD), 0.43 (CZ) and 0.37 ('flat'), respectively. 
Using the ERBL evolution equations we can compute the second moments at the scale 
/x 2 = 4 GeV 2 . We find the values 0.20 (asymptotic), 0.24 (AdS/QCD) and 0.38 (CZ). 
The agreement between the AdS model value and the lattice results is better than the 
result found for the asymptotic model and CZ model. We also note that the measurement 
of the pion DA in diffract ive di-jet production reported by the E791 Collaboration [58] 
supports a centrally-peaked DA such as the asymptotic and AdS/QCD models. The 
second moment alone, while providing important information for the pion DA, will not 
put strong constraints on the shape of the pion DA, since it is a quantity obtained by 
integrating the DA over the whole range of x. 



III. PION-PHOTON TRANSITION FORM FACTORS 

The pion-photon transition form factor can be extracted from the two-photon process 
7*(?i)7*(?2) —> tt°- When both photons are off-shell with virtuality Q\ = —q\ and 
Q\ = — gf, the form factor is denoted as F 7T1 *(Q\ ) Q|)- ^ n ^he case °f one photon being 
on mass-shell the form factor is denoted as F^^Q 2 ). 



A. Leading order results 



Brodsky and Lepage [71 HU] predicted the behavior of F 7r7 (Q 2 ) at leading order of 
a s (Q 2 ) and leading twist as 



V o Jo x 



rn 



l + 0(a s , Q2 



(23) 



where x is the longitudinal momentum fraction of the quark struck by the virtual photon 
in the hard scattering process and x = 1 — x is the longitudinal momentum fraction of 
the spectator quark. It was argued [7j that the boundary condition of DAs vanishing at 
the end-points faster than x e for some e > would enable one to replace 0(x, xQ) by 
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(x,Q) since the difference is non-leading 2 , and Eq. (23) becomes 



X 



1 + [a s , 



m 



(24) 



The replacement is sound when one is interested in the leading order behavior of the 
TFF and particularly for the behavior at the asymptotic limit Q 2 — >■ oc, which is one 
of the main purposes of Ref. [7J. However, this approximation is not justified for the 
calculation at finite Q 2 region where one needs to take into account the evolution effects 



and NLO corrections. The dominant contributions to the integrals in Eqs. (23) and 



(24) come from small x region, e.g., 3/4 of the contributions coming from x < 0.5 for 



the asymptotic DA. At the same time the evolution changes the shape of the DA more 



significantly in the small x region. Thus the calculations with Eq. (23) involve a much 



less evolved DA than Eq. (24). The difference between the calculations using Eqs. (23) 



and (24) could be sizable. Unfortunately, Eq. (24) has been widely used in the literature 



as the starting point to calculate high-order corrections to the pion TFF, see e.g., [201424] . 
Similar replacement has been done in the study for other exclusive processes. 

It is essential to consider the transverse momentum dependence in both the hard- 
scattering amplitude and the LFWF in order to describe the data at finite Q 2 [T6] . 
Taking into account the fc^-dependence, the pion-photon transition form factor is given 
by [HITS] 



dx 



oo j2 



d 2 k ± 
16tt 3 



T H (x, Q 2 , k±)ip qil/n (x, k|_), 



(25) 



where 



T H (x,Q 2 ,k ± ] 



q± • (xqx + k ± ) 

— + [x x] , 



(26) 



cfi_(xqx + k ± f 

is the hard scattering amplitude and = Q 2 . Using Eq. (25) and a Gaussian type 



LFWF one can reproduce the curve displayed by the experimental data at low Q 2 [IB] . 
With Eqs. (23) and (24) the calculations will be near constant for all Q 2 . 



It was actually pointed out that the replacement xQ — > Q/2 is more appropriate. 
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Musatov and Radyushkin have shown [59] that if LFWF depends on the transverse 
momentum only through k^, Eq. (25) can be simplified as 



4 r dr f x ^ d 2 k i 



(27) 



For the model wavefunction Eq. ^ we can factor out the Q 2 dependence of the DA at 



low Q 2 , Eq. (11), and include the QCD evolution for higher momenta through the ERBL 
solution of the DA. One obtains 3 [59] 



4 



V o Jo x 



1 — exp 



xQ 2 
2n 2 x 



(28) 



The pion TFF depends on Q 2 through the exponential factor and the pion DA. Since 
we have explicitly factored out the low (3 2 -dependence, the distribution amplitude in 



Eq. (28) contains only the hard ERBL evolution. The exponential factor is important, 



especially for small x and small Q 2 , thus it controls the curvature of Q 2 F 7T1 (Q 2 ) vs. Q 2 
at low Q 2 . The behavior of the pion TFF at high Q 2 is determined dominantly by the 
pion 'hard' DA which should evolve in a logarithmic manner. The exponential factor also 
plays a role to regularize the calculation with the 'flat' DA, which otherwise involves a 
divergent integral. 



Inserting Eq. (16) into Eq. (28) we can write the transition form factor as 



oo 



n=0,2,4- 



a n (po) 



[ &xxCl /2 {2x- 1) 

JO 





7n/A) " 


Uw)J 





1 — exp 



xQ 2 
2k 2 x 



(29) 



which displays the soft and hard dependence. We need to set xQ = /x for xQ < /x , 



which assures the convergence of Eq. (29). Equations ([28]) and (29) clearly show that 
the pion TFF at any given Q 2 is determined by /x ) and all evolved DAs from /x 
to Q, with the less-evolved DAs providing major contributions. For example, half of the 



3 Enforcing the cut-off < x(l — x)Q 2 for the soft LFWF discussed in Section I, Eq. (28) becomes 

X - ex P ("2&) 



where Q = xxQ with x = mm(x,x). The two 



expressions coincide for x < 0.5 and the differences are negligible unless Q 2 < 1 GeV 2 . However, for 
other exclusive processes that are sensitive to the large-x region the difference may be sizable. 
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FIG. 5: The pion-photon transition form factor shown as Q 2 F 7T7 (Q 2 ) calculated using Eq. (28) 
with different prescriptions for solid curve - AdS (x, xQ), dashed curve - AdS (x), 

thick-dashed curves - AdS (x,(3), and dash-dotted curve - asy (x). P qq = 0.50. The data are 
taken from [TJ |5l E] . 

contributions at Q — 3 /x come from /x ) for the asymptotic DA, and this ratio is 
much higher for broad models for the pion DA. The contributions from 0(x,/x o ) remain 
significant even when Q ~ 5/x . Thus the evolution effect hardly shows up until Q 2 is 



very large. On the other hand, if one uses the distribution Q) in Eqs. (28) and (29) 
the evolution effect will be overestimated. 



We compare results calculated using 0(x,/x o ), (j)(x,xQ) and Q) in Eq. (28). The 
results for the AdS and CZ models for the pion DA are shown in Figs. [5] and [6j respectively. 
The valence probability P q q = 0.50 has been adopted in these calculations. Using Q) 
will unjustifiably reduce the predictions substantially for Q 2 > 10 GeV 2 . We conclude 
that the evolution effect at leading order will not bring any large corrections to the 
calculation for the pion transition form factors. It is a good approximation to use /i ) 
in the pQCD calculation for exclusive processes. This conclusion can be expected to hold 
when the evolution is considered at NLO as well. 



Analytical expression exists for each term in the sum in Eq. (29), though the expression 
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FIG. 6: Similar as in Fig. [5] but for the CZ model for the pion DA. 

becomes extraordinary long and tedious for large n. The first term corresponds to the 
results with the asymptotic DA, 



Q*Fi?m = 2/„^_ ^ _ g e |&r[0, g] j , (30) 
where T[0,x] is the incomplete gamma function. At the asymptotic limit Q 2 oo, Eq. 



>>2 T?AS(r\2\ n r Q 



(30) gives Q F^(Q ) —> 2f n as expected. A slightly more complicated expression exists 



for the CZ model for the pion DA. 



B. Next-to- leading Order Corrections 

The next-to-leading order corrections have been studied [201424] under the assump- 
tion of cj){x,xQ) ~ Q), using the standard hard scattering approach when the k±- 
dependence in the hard scattering amplitude is ignored. As illustrated in the last section, 
a properly treatment of evolution is required. So it is necessary to revisit the NLO cal- 
culations with cj>{x, xQ). Assuming that the fc^-dependence of the LFWF introduces the 
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FIG. 7: Effect of NLO corrections on the pion to photon transition form factor Q 2 F 7T7 (Q 



The curves without markers are the results calculated using Eq. (28) with <j)(x) for the four 



models of the pion DA: solid curve - CZ model, dashed curve - 'flat' model, thick-solid curve - 
AdS model, and dash-dotted curve - asymptotic model. The curves with markers are the NLO 



results calculated using Eqs. (31) and (32). P q q = 0.50. The data are taken from [U|5l|6j. 



same exponential factor for the small Q region 4 the TFF can be expressed as 



where 



Q 2 F^°(Q 2 ) 



V3 



dx Th(x : Q 2 ) (j)(x : xQ) 



1 — exp 



xQ 2 
2n 2 x 



2x = 1 , OLs^r)^ 1 
X 47T X 



T H (x,Q' z ) = ^ + ^f^C F 



—9 \nx + ln 2 x 

x 

+ (3 + 21nx)ln ( % 



(31) 



(32) 



and (j){x,xQ) is the 'hard' DA evolved at the next-to-leading order [60J, except for the 
'flat' DA which cannot be evolved as discussed in Section 2. The regularization scale is 
commonly taken as Hr = Q to eliminate otherwise large logarithm terms. 



Dropping this exponential factor will hardly affect the calculation at large Q 2 . 
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The numerical results for the pion TFF with the four DA models discussed in Section 
II are shown in Fig. [7| The NLO corrections vary according to the models used for the 
pion DA. The corrections at Q 2 ~ 5 GeV 2 are about 20%, 17%, 11%, and 15% for the 
asymptotic, AdS, CZ, and 'flat' models. The corrections at Q 2 ~ 40 GeV 2 are still more 
than 7% for all the DAs. Thus it is necessary to take into account these corrections for 
a large range of Q 2 . 

C. Higher Order and Higher Fock State Contributions and Dependence on P q q 

The calculations for the transition form factors depend on the non-perturbative input, 
i.e. the soft LFWF. Consequently, it has been long argued that the pion-photon transition 
form factor is a particularly suitable process in determining the pion LFWF and DA. 
As discussed in Section II, we have treated the only parameter k in the model LFWF 
(Eq. ([6])) as a parameter which is constrained, though not very strictly, by the probability 
of finding the valence Fock state and the mean square transverse momentum. In the above 
calculations we have adopted P q q = 0.50. The next-to-leading order predictions with 
P q q = 0.50 for the pion-photon transition form factor are smaller that the experimental 
data, particularly for the Q 2 < 10 GeV 2 region. To improve the agreement between the 
calculations and experimental data one could use a larger value for P q q. For example, 
using Pqq = 1.0 will give a much better agreement for the calculations with the CZ model 
for the pion DA. However, a much larger value of P q q than 0.5 will result in a much 
smaller value for the root of mean square transverse momentum of the valence quarks 
compared to the value obtained from the charge radius of the pion. 

It is shown in [24] that the next-to-next-to-leading corrections are much smaller that 
the next-to-leading corrections. However, the contributions from higher Fock states (e.g., 
\qqqq)) are important at low Q 2 . Figure [8] (b) illustrates such a contribution where 
each photon couples directly to a qq pair. Such higher-twist contributions $^ e .^ e 
are necessary to derive the low energy amplitude for Compton scattering — >> ^H, 
which is proportional to the total charge squared e 2 H = (e* + ej) 2 of the target. These 
contributions are suppressed by the factor [1/Q 2 ) n at large Q 2 , where n can be understood 
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FIG. 8: Leading-twist contribution (a) and a possible higher-twist contribution (b) to the 



process 77* tv°. 



as the number of qq pairs in the higher Fock states. An analysis of these contributions 
using the framework of AdS/QCD is presented in [50J. To estimate these higher Fock 
state contributions we adopted a phenomenological model as in [27] 

^ 7 (0)/2 



(1 + Q 2 /A 2 ) 2 ' 



(33) 



where F^O) = l/(47r 2 / 7F ) is the PCAC result and A can be treated as a parameter. The 
contributions are less than 1% for Q 2 > 10 GeV 2 and thus can be safely ignored. 

The total contribution from the valence Fock state and the higher Fock states is the 
sum of Eqs. (31) and (33). The results calculated with the choice of P q q = 0.5 and 



A = 1.1 GeV in Eq. (33) are compared with the data in Fig. [9j The agreement at the 
low Q 2 region is vastly improved due to the inclusion of higher Fock state contributions. 
However, the higher Fock state contributions are negligible for Q 2 > 10 GeV 2 and it 
is in this large Q 2 region that the four models of the pion DA, discussed in Section II, 
give very different predictions for the Q 2 dependence of the pion-photon transition form 
factor. The results with the asymptotic DA are smaller than the BABAR data and, as 
expected, do not exhibit a strong Q 2 dependence. The results with the 'flat' DA show 
a substantial and continuous growth with Q 2 , which is in disagreement with the QCD 
prediction that the pion TFF should approach its asymptotic value of 2f n at Q 2 — >■ 00. 
In fact, one cannot apply the ERBL evolution equation to the 'flat' DA since it does not 
satisfy the boundary condition of the pion DA vanishing at x = and x — 1. The results 
with the AdS model and CZ model for the pion DA lie in between the predictions of 
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FIG. 9: The tt° — 7 transition form factor including contributions from the valence Fock state 
(Eq. ([31])) and higher Fock states (Eq. ([33])) of the pion. The data are taken from [HO [6]. 



the asymptotic DA and the 'flat' DA. The results with the CZ DA show a fast growth 
with Q 2 compared with the AdS DA over the range of 10 GeV 2 < Q 2 < 100 GeV 2 . 
We note that the BABAR data for Q 2 > 20 GeV 2 suffer larger uncertainties as compared 
with the low- and medium-Q 2 regions. We also note that the 'flat' and CZ models of 
the pion DA will produce much larger values for the /y-photon and //-photon transition 
form factors than the results reported by the BABAR Collaboration for Q 2 > 15 GeV 2 
[31 [4] and at Q 2 = 112 GeV 2 [61J. Figure [9] also shows that the calculations approach 
the asymptotic limit value Q 2 F 7T7 (Q 2 — >• 00) = 2f n very slowly since the DA evolution 
introduces a logarithm (3 2 -dependence via [In (Q 2 /Aqcd) /WMoMqcd)]" 7 ™- 

We investigate the dependence of our calculations on P qq by allowing P qq to be in the 



range of 0.5 ~ 0.8. The valence Fock state contributions calculated with Eq. (31) are 
One can see that the calculations for Q 2 > 30 GeV 2 depend on P c 



shown in Fig. 



10 



qq 



very weakly, though the dependence at the lower Q 2 region is much more significant. The 
four models of the pion DA give very different predictions for the pion-photon transition 



for the region of Q 2 > 30 GeV 2 , regardless the value of the P c 



qq- 
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FIG. 10: The tv° — 7 transition form factor calculated with Eq. (31) for P qq = 0.5 and P qq — 0. 
The data are taken from [H O [6]. 



It is very difficult to accommodate the BABAR large-Q 2 data with the QCD calculations 
using the asymptotic, AdS, and CZ models for the pion DA. The calculations with the 
'flat' model of the pion DA can produce a rapid growth for the pion TFF shown by 
the BABAR data. However, the calculations with the same DA model underestimate 
significantly the pion TFF at the low Q 2 ) and the prediction for the pion TFF at the 
asymptotic limit Q 2 — » 00 violates seriously the Brodsky-Lepage limit of Q 2 F 7T7 (Q 2 
00) = 2f v . 

It was pointed out in [38] that using a contact interaction for the quark-antiquark in- 
teraction in the Dyson- Schwinger equations (i.e. treating the pion as a point-like bound 
state) produces a 'flat' DA and gives predictions for the pion electromagnetic form fac- 
tors [62J and transition form factor that are in striking disagreement with completed 
experiments. In Ref. [3]] the pion is treated as an elementary field in the triangle graph 
and the simple expression obtained as F 7T o 1 (Q 2 ) ~ ^(ln^) 2 (with m = 132 MeV) is 
able to reproduce the BABAR data for the pion TFF. We would like to emphasize that 
although the chiral field theory is a useful approximation for some long-wavelength, soft 
processes, it is inapplicable to the hard scattering regime of the BABAR data. In fact, 
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the compositeness of the pion in terms of quarks and gluons has been verified in high 
energy experiments both in inclusive reactions (such as the Drell-Yan process for pion- 
nucleon collisions) and many hard exclusive reactions (such as the pion form factor at 
large spacelike and timelike momentum transfers and large angle scattering processes 
such as 77 — >■ tttt and pion photoproduction) . It is also not necessary to treat the pion as 
elementary to prove chiral anomalies or the Gell-Mann-Oakes-Renner (GMOR) relation. 
Such relations are standard consequences of QCD for a composite pion [63]. Employing 
the BMS models [25] for the pion DA, which were determined utilizing the CELLO and 
CLEO data for the pion TFF, will produce a (^-dependence for the pion TFF similar to 
that obtained with the AdS model for the pion DA. A recent analysis [39] of all existing 
data (CELLO, CLEO and BABAR ) performed in a framework similar to [25] suggest that 
it is not possible to accommodate the high-Q 2 tail of the BABAR data with the same 
accuracy as the analysis of the CELLO and CLEO data. 

We note that there are several theoretical studies [271436] trying to reproduce the 
BABAR data for the pion TFF, apart from those using the 'flat' form for the pion DA 
[T3l - fT5l [26] . It was claimed in [271 EE] that a much broader DA than the asymptotic 
form (but which still vanishes at the end-points) would be able to explain the BABAR 
results. The Regge approach was employed in [29, 30j to explain the BABAR data. On 
the other hand, there are also theoretical calculations suggesting that the BABAR data 
are not compatible with QCD calculations [37-40J. 

We would like to remark that more accurate measurements of the pion-photon tran- 
sition form factor at the large Q 2 region will be able to distinguish the various models of 
the pion DA under discussion. 

D. The transition form factor for the pion- virtual-photon 

The above analysis can be easily extended to the case in which the photons involved 
are both off mass-shell, i.e., for the form factor F 7T1 *(Q\ ) Q\) ) by replacing the hard- 
scattering amplitude T H with the corresponding expression. At leading order T H has the 



24 



form 



T-^-\ X , Ql M = ^_, ,34) 



where x = 1 — x. For the expression at next-to-leading order we refer the readers to 
reference [22] . 

A significant difference from the case with a real and a virtual photons where X# = 



l/(xQ 2 ) is that Eq. (34) is not divergent at the end-points. Thus considering the /in- 
dependence will not bring as large corrections as for F 7T1 {Q 2 ) ) and the transition form 
factor F^* is much less sensitive to the end-point behavior of the pion DA than F^. 

We note that the kinematic region satisfying Q\ = Q\ is particularly interesting since 
in this region the amplitude Tjj becomes independent of x and thereby the transition 
form factor is largely described by the normalization of the pion DA, which is model and 
Q 2 independent. Ignoring the weak Q 2 dependence introduced by the consideration of 
k± dependence and the NLO corrections in a S) which are both expected to be small at 
large Q 2 , we have 

QlF^AQl Ql) -+ \U for Qi = Q 2 > a few GeV. (35) 



We make the remark that Eq. (35) is expected to work for the range Q\ = Q\ ~ 10-20 
GeV 2 which is accessible by the current experiments. So measurements of F^{Q\^Q\) 
under these conditions would provide another test of pQCD analysis of exclusive pro- 
cesses. 



The numerical results for F 7T1 ^{Q\ ) Q\) calculated at NLO are given in Fig. 11 for 
Q\ = 2 GeV 2 . The four models give similar predictions for Ql up to 15 GeV 2 . At large 
Ql the results with the CZ model are much larger compared with the asymptotic and 
AdS models. We found that the NLO corrections at this range of Q 2 are less than 10% 
for the asymptotic and AdS models while the corrections to the CZ and 'flat' models are 
negligible. The higher-twist effects at this range of Q 2 could be expected to be minimal. 
Thus the difference on the prediction for this transition form factor is a direct reflection 
of different behavior of the pion DA. Measurements of this form factor at the kinematic 
region Ql ~ 20 Q\ with Q\ being about a few GeV 2 would provide a good laboratory to 
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FIG. 11: The doubly virtual transition form factor for Q\ — 2 GeV 2 , calculated with Eqs. (31) 



and (34). The solid, dashed, thick-solid, and long-dashed curves are the results with the CZ, 



'flat', AdS, and asymptotic models for the pion DA. 



distinguish the middle-peak DA, such as the asymptotic and AdS models, from the CZ 
model. 



IV. THE 77-PHOTON AND //-PHOTON TRANSITION FORM FACTORS 

According to the SU(3)i? quark model, the three charge neutral states in the nonet of 
pseudoscalar mesons are 7r°, rjs and 771 . The latter two mix to give the physical particles 77 
and 7/. It can be expected that the states 7r°, ?y 8 and 771 have the same form of distribution 
amplitude (and the same fc^-dependence in the light-front wavefunctions) , 

<t>pix) = fpfi?), (36) 

with P denoting 7T°, rjs and 771, and fp being the corresponding decay constant. 
The transition from factors for the 7r°, rjs and rji can be expressed as 

4 f 1 

Q 2 F Pl {Q 2 ) = ^c P dxT H (x,Q 2 )0 P (x,xQ) 
v3 Jo 
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1 — exp 



xQ 2 
2n 2 x 



(37) 



where cp = 1, and ^2 for 7r°, 7y 8 and 771, respectively and T H {x ) Q 2 ) is given by 



Eq. (32) at next-leading order of QCD running coupling constant. 



The transition form factors for the r] and 7/ result from the mixing of F mi and F mi , 

(38) 



cos 9 — sin 9 


11 


F \ 


sin 9 cos 9 j 







where 9 is the mixing angle which has been the subject of extensive studies [64J. In this 
work we adopt 9 = -14.5° ± 2°, f 8 = (0.94 ± 0.07)/^ and f ± = (1.17±)f n [65J. The 
same value of n has been used for the three charge neutral states and A = 1.1 GeV is 



adopted in Eq. (33). The results for the 77-photon and //-photon transitions form factors 



are shown in Figs. [12] and [13] respectively. The data favor the AdS and the asymptotic 
models for the meson DA. One may fine-tune the parameters /§, and fi to make the 
calculations with the asymptotic form and the AdS form to give better agreement with 
the data. However it is almost impossible to make the calculations with the CZ form and 
the 'flat' form to describe the data in both the low- and high-Q 2 regions simultaneously 
for the two transition form factors, although these two forms are favored to explain the 
rapid growth of the BABAR data [lj for the pion-photon transition form factor at large 

Q 2 - 

The DA models that could explain the BABAR measurements for the pion-photon 
transition form factors at large values of Q 2 will fail in QCD calculations for the other 
processes, including the 77-photon and //-photon transition form factors reported by the 
BABAR Collaboration. This may suggest that the BABAR measurements at large Q 2 are 
not a true accurate representation of the pion-photon transition form factor, a perspec- 
tive that has been suggested in [371 EEj. This may also indicate that there are some 
inconsistencies among the results for the transition form factors of the 7T°, tj and 7/ 

In a recent paper [40J, Wu and Huang have studied the dependence of the photon-to- 
meson transition form factors on the model parameters (including quark masses, mixing 
angle, as well as an intrinsic charm component) in a light-front perturbative approach. 
It is found that the agreement of the predictions of their model with the experimental 
data can be somewhat improved by adjusting these parameters within their reasonable 
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FIG. 12: The 77 — 7 transition form factor shown as Q 2 F r]1 (Q2). The thick-dashed, thick-solid, 
thin-dashed, thin-solid curves are the results calculated with the asymptotic, AdS, CZ and 'flat' 
models for the meson DAs, respectively. Data are taken from [6J (CLEO) and [Sill] (BABAR). 
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FIG. 13: The rf — 7 transition form factor shown as Q 2 F r] / 1 (Q2)- The thick-dashed, thick-solid, 
thin-dashed, thin-solid curves are the results calculated with the asymptotic, AdS, CZ and 'flat' 
models for the meson DAs, respectively. Data are taken from [6J (CLEO) and [21 El (BABAR). 
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regimes, but the data for the pion-photon transition form factor in the entire Q 2 domain 
cannot be explained consistently. 

Recently Kroll [28] analyzed the n — 7, 77 — 7, and 7/ — 7 transition form factors using 
the modified hard scattering approach in which the transverse-momentum-factorization 
is combined with a Sudakov factor. The distribution amplitude of mesons is constrained 
to contain the first three nontrivial terms in the Gegenbauer expansion for the DA, and 
a Gaussian form is assumed for the fc^-dependence in the wave function. By adjusting 
the three parameters - the two coefficients, a 2 and in the Gegenbauer expansion, and 
the transverse size parameter op, reasonably good agreement with experimental data 
was achieved. The best fit for the pion-photon transition form factor presented in Fig. 
4 of [28J is very similar to our results calculated with the AdS model for the pion DA 
(see Fig. [9]). We note that in order to describe the three transition form factors, very 
different values for the three parameters are chosen in Ref. [28j for the 7r , r]s an d 771 : 
al = 0.20, a\ = 0.01, a* = 0.40 GeV" 1 ; a 8 = -0.06, a\ = 0, a 8 = 0.84 GeV" 1 ; 
a\ = —0.07, a\ = 0, and a 1 = 0.74 GeV -1 . Such a choice of parameters suggests a very 
large SU(3)i? symmetry breaking between the DAs of the tt° and rf and a very little 
SU(3)i? symmetry breaking between the rjs and 771 . Our view is that this remains as a 
possibility, but it is unlikely since the ttq and rjs belong to the octet and 771 belongs to the 
singlet of the pseudoscalar mesons. Furthermore, there is no evidence supporting such 
a large SU(3).f breaking in other processes, e.g., decay processes involving pseudoscalar 
mesons [64J. In fact, the CLEO's measurements of the meson-photon transition form 
factors [6] suggested that the tt° and rjs have very similar non-perturbative dynamics and 
thereby similar light- front wavefunctions and distribution amplitudes. 



V. SUMMARY 



The photon-to-meson transition form factor measured in 7*7 — »> M is the simplest 
hadronic amplitude predicted by QCD. Measurements from electron-positron colliders 
provide important constraints on the non-perturbative hadron distribution amplitude, 
a fundamental gauge-invariant measure of hadron structure. The meson distribution 
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amplitude Q) evolves in pQCD according to the ERBL evolution equation, which 
is based on first-principle properties of QCD. Important constraints on the distribution 
amplitude have been obtained using lattice gauge theory. We have analyzed in detail 
four models for the 7r°, 77, and rf distribution amplitudes that have been suggested in the 
literature, including their QCD evolution with logQ 2 . 

We have calculated the meson-photon transition form factors for the 7T°, 77 and 7/, 
taking into account effects which are important for the calculations at finite Q 2 . These 
effects include the fc^-dependence of the hard-scattering amplitude and light-front wave- 
functions, the evolution effects of the pion distribution amplitude, and NLO corrections 
in a s . We have pointed out that a widely-used approximation of replacing <j){x,xQ) with 
<f){x,Q) in the hard-scattering formalism will significantly, and unjustifiably, reduce the 
predictions for the magnitude of hard exclusive amplitudes. 

It is found that in order to explain the experimental data at Q 2 < 10 GeV 2 one needs 
to take into account the contributions from higher Fock states of the mesons, although 
these contribution are negligible for the larger Q 2 region. The four models of the meson 
DA discussed in this article give very different predictions for the Q 2 dependence of 
the meson-photon transition form factors in the large Q 2 region. The predictions based 
on the AdS/QCD and light-front holography for the pion distribution amplitude agree 
well with the experimental data for the 77- and //-photon transition form factors, but 
they disagree with the data for the pion-photon transition form factor reported by the 
BABAR Collaboration. The calculations with the CZ model agree with the BABAR data 
for the pion-photon transition form factor reasonably well, but the predictions are much 
larger than the data from the CLEO and BABAR Collaborations for the rj- and //-photon 
transition form factors. The calculations with the 'flat' distribution amplitude, which has 
been advocated in explaining the BABAR large-Q 2 data for the pion transition form factor, 
disagree strongly with the CLEO and BABAR data for the 77- and //-photon transition form 
factors. 

We investigated the dependence of the calculations on the probability of valence Fock 
state of the pion P q q. It was found that the four models of the meson DA give very 
different predictions for the meson-photon transition form factor in the region of Q 2 > 30 
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GeV 2 for P q q to be in the reasonable range of 0.5 ~ 0.8. More accurate measurements of 
the meson-photon transition form factor in the large Q 2 region will be able to distinguish 
the four commonly used models of the pion DA. 

The BABAR data for the pion-photon transition from factor exhibit a rapid growth at 
high Q 2 , but this feature is missing for the 77- and ?/-photon transition form factors. The 
rapid growth of the large-Q 2 data for the pion-photon transition form factor reported by 
the BABAR Collaboration is difficult to explain within the current framework of QCD. 
This is a viewpoint first expressed by Roberts et al. [38] in their Bethe-Salpeter/Dyson- 
Schwinger analysis of the pion-photon transition form factors. If the BABAR data for the 
meson-photon transition form factor for the 7T° is confirmed, it could indicate physics 
beyond-the-standard model, such as a weakly-coupled elementary C = + axial vector or 
pseudoscalar z° in the few GeV domain, an elementary field which would provide the 
coupling 7*7 —> z° —> 7T° at leading twist. We would like to remark that a high-mass 
state of about 10 GeV has been envisaged in [29j to explain the BABAR data for the 
pion TFF [66]. We thus emphasize the importance of additional measurements of the 
meson-photon transition form factors. 
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